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The electric field and the hole and electron concentrations are found for 
reverse biased junctions in which one side is either intrinsic (I) or so wealdy 
doped that the space charge of the carriers cannot he neglected. The analysis 
takes account of spare charge, drift, diffusion and non linear recombination. 
A number of figures illustrate the penetration of the electric field into a PIN 
structure with increasing bias for various lengths of the I region. For the 
junction between a highly doped and a weakly doped region, the reverse cur- 
rent iiicreases as the square root of the voltage at high voltages; and the space 
charge in the weakly doped region approaches a constant value that depends 
on the fixed charge and the intrinsic carrier concentration. 

The mathematics is greatly simplified by expressing the equations in 
terms of the electric field and the sum of the hole and electron densities. 

I. INTRODUCTION 

Applications have been suggested for semiconductor structures having 
both extrinsic and inti'insic regions. Examples are the "swept intrinsic" 
structure, in which a region of high resistivity is set up by an electric 
field that sweeps out the mobile carriers, and the analogue transistors, 
where the intrinsic region is analogous to the vacuum in a vacuum tube. 
However, the j unction between an intrinsic region and an iV^ or P region 
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is considerably less well understood than the simple NP junction. Most 
of the assumptions that make the NP case relatively simple to deal with 
do not apply to junctions where one side is intrinsic. Specifically, the 
space charge is that of the mobile carriers ; thus the flow and electrostatic 
problems cannot be separated as they can in PA'' junction under reverse 
bias. The following sections analyze the A'"-intrinsie - P structure under 
reverse bias. 

For a given material with fairly highly doped extrinsic regions, the 
problem is defined by the length of the intrinsic region and the applied 
voltage. Taking the intrinsic region infinitely long gives the solution for 
a simple A''-intrinsic or P-intrinsic structure. The results are given and 
plotted in terms of the electric field distribution. From this the potential, 
space charge and carrier concentrations can be found; so also can the 
current-voltage curve. The final section considers the case where the 
middle layer contains some fixed charge but where the carrier charge 
cannot be neglected. 

Qualitative Discussion of an N-intrinsic-P Structure 

Consider an iV-intrinsic-P structure where the intrinsic, or /, region is 
considerably wider than the zero bias, or built-in, space charge regions 
at the junctions, so that there is normal intrinsic material between the 
junctions. The field distribution at zero bias can be found exactly from 
the zero-current analysis of Prim.' Throughout the intrinsic region, hole 
and electron pairs are always being therinally generated and reeombining 
at a rate determined by the density and properties of the traps, or recom- 
bination centers. Under zero bias the rates of generation and recombina- 
tion are everywhere equal. Suppose now a reverse bias is applied causing 
holes to flow to the right and electrons to the left. Some of the carriers 
generated in the intrinsic region will be swept out before reeombining. 
This depletes the carrier concentration in the intrinsic region and hence 
raises the resistivity. It also produces a space charge extending into the 
intrinsic layer. The electrons are displaced to the left and the holes, to 
the right. Thus the space charge opposes the penetration of the field 
into the intrinsic region; that is, the negative charge of the electrons on 
the left and positive charge of the holes on the right gives a field 
distribution with a minimum somewhere in the interior of the intrinsic 
region and maxima at the NI and IP junctions. If holes and electrons 
had equal mobilities, the field distribution would be symmetrical with a 
minimum in the center of the intrinsic region. Likewise, the total carrier 

iR. C. Prim, B. S. T. J., 32, p. 665, May, 1953. 
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concentration (holes plus electrons) would be symmetrical with a maxi- 
mum in the center. As the applied bias is increased the hole and electron 
distributions are further displaced relative to one another and the space 
charge increases. Finally, at high enough biases, so many of the carriers 
arc swept out immediately after being generated that few carriers are 
left in the intrinsic region. Now the space charge decreases with increas- 
ing bias until there is negligible space charge, and a relatively large and 
constant electric field extends through the intrinsic region from junction 
to junction. This may happen at biases that are still much too low to 
appreciably affect the high fields right at the junction or in the extrinsic 
layers, which remain approximately as they were for zero bias. 

The current will increase with voltage until the total number of 
carriers in the intrinsic region becomes small compared to its normal 
value. After that, there is negligible further increase of current with 
voltage. All the carriers generated in the intrinsic region are being swept 
out before recombining. In general, the current will saturate while the 
minimum field in the intrinsic region is still small compared to the 
average field. 

Comparison with the NP Structure 

The analysis is more difficult than in a simple reverse-biased NP 
structure. In the NP case there is a well defined space charge region in 
which carrier concentration is negligible compared to the fi-ted charge of 
the chemical impurities; so the field and potential distributions are easily 
found from the known distribution of fixed charge. Outside of the space 
charge region are the diffusion regions in which the minority carrier con- 
centration rises from a low value at the edge of the space charge region 
to its normal value deep in the extrinsic region. However, there is no 
space charge in this region because the majority carrier concentration, 
by a very small percentage variation, can compensate for the large per- 
centage variation in minority carrier density. The minority carriers flow 
by diffusion. Since the disturbance in carrier density is small compared to 
the majority density, the recombination follows a simple Unear law 
(being directly proportional to the excess of minority carriers). Thus 
the minority carrier distribution is found by solving the simple diffusion 
equation with linear recombination. 

None of these simplifications extend to the NIP or NT or IP structure. 
There is, in the intrinsic region, no fixed charge; hence the space charge 
is that of the carriers. There is no majority carrier concentration to 
maintain electrical neutrality outside of a limited space charge region. 
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It is necessary to take account of (1) space charge, (2) carrier drift, (3) 
carrier diffusion and (4) recombination according to a nonlinear bi- 
molecular law. Of these four, only space charge and recombination are 
never simultaneously important in practical cases. Nevertheless certain 
simplifications can be made if the problem is formulated so as to take 
advantage of them. The field and carrier distributions in the intrinsic 
region are found by joining two solutions: one solution is for charge 
neutrality; the other, which we shall call the no-recombination solution 
is for the case where the recombination rate is negligible compared to the 
rate of thermal generation of hole electron pairs. We shall show that in 
practical cases the ranges of validity of the two solutions overlap; that 
is, wherever recombination is important, we have charge neutrahty. 

Prim's Zero-Current Approximation 

Prim* derived the field distribution in a reverse biased NIP structure 
on the assumption that the hole and electron currents are negligibly small 
differences between their drift and diffusion terms, as in the zero-bias 
case. He showed that the average diffusion current is large compared 
to the average current. However, as it turns out, this is misleading. 
Throughout almost all of the intrinsic region (where the voltage drop 
occurs in practical cases) the diffusion current is comparable to or 
smaller than the total current. The larger average diffusion current comes 
from the extremely large diffusion current in the small regions of high 
space charge at the junctions. Prim's analysis, in effect, neglects the space 
charge of the carriers generated in the intrinsic region. These may be 
neglected in calculating the field distribution if the intrinsic region is 
sufficiently narrow or the reverse bias sufficiently high. In the appendix 
we derive the limits within which Prim's calculation of the field and 
potential will be valid. The range will increase with both the Debye 
length and the diffusion length in the intrinsic material. However, in 
cases of practical interest the zero-current approximation may lead to 
serious errors in the field distribution and give a misleading idea of the 
penetration of the field into the intrinsic region. The present, more 
general analysis, reduces to Prim's near the junctions where the zero- 
current assumption remains vaUd. The zero current approximation was, 
of course, not intended to give the hole and electron distributions in the 
intrinsic region or to evaluate the effects of interacting drift, diffusion 
and recombination. 

*Ibid. 



THEORY OF THE SWEPT INTRINSIC STRUCTURE 1243 

Outline of the Following SecHons 

Sections II through V deal mth the ideal case of equal hole and elec- 
tron mobilities. Here the problem is somewhat simplified and the physics 
easier to visualize because of the resulting sj'mmetry. In Section VI, 
the general case of arbitrary mobilities is solved by an extension of the 
methods developed for solving the ideal case. The technique is to deal 
not with the hole and electron flow densities but with two linear com- 
binations of hole and electron flow densities that have a simple form. 

Section II deals with the basic relations and in particular the fonnula 
for recombination in an intrinsic region for large disturbances in carrier 
density. The nature and range of validity of the various approximations 
are discussed. Section III derives the field distribution in regions where 
recombination is small compared to pair generation. Section IV treats 
the recombination region and the smooth joining of the recombination 
and no-recombination solutions. Section V considers the role of diffusion 
in current flow and the situation at the junctions where the field and 
carrier concentration abruptlj'^ become large. The change m form of 
the solution near the junctions is shown to be represented by a basic in- 
stability in the governing differential equation. Section VI extends the 
results to the general case of unequal mobilities. Section VII deals with 
the still more general ease where there is some fixed charge in the "in- 
trinsic" region. If the density of excess chemical impurities is small com- 
pared to the intrinsic carrier density, the solution remains unchanged in 
the range where recombination is important. In the no- recombination 
region the solution is given by a simple first order difi"erentiatial equation 
which can be solved in closed form in the range ^^■here the carrier flow is 
by drift. The fixed charge may have a dominant effect on the space 
charge even when the excess density of chemical impurities is small com- 
pared to the density n,- of electrons in intrinsic material. Consider, for 
example, a junction between an extrinsic P region and a weakly doped 
n region having an excess density N = Nd — Na of donors. In the limit, 
as the reverse bias is increased and the space charge penetrates many 
diffusion lengths into the ?i region, the field distribution becomes linear, 
corresponding to a constant charge density equal to 

h[N + y/N^ + 8 n.-=£VL/] 

where L,- is the diffusion length in the weakly doped n type region and £. 
is the Debye length for intrinsic material. For germanium at room tem- 
perature £/Li is the order of 10~ . Thus, in this example, a donor density 
as low as lO" cm^^ will have an appreciable effect on the space charge. 
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II. BASIC RELATIONS 

The problem can be stated in terms of the hole density p, the electron 
density n, and the electric field E and their derivatives. Let the distance 
X be measured in the direction from A'' to P. The field will be taken as 
positive when a hole tends to drift in the -\-x direction. The field in- 
creases in going in the +x direction when the space charge is positive. 
Poisson's equation for intrmsic material is 

^ = aip- n) (2.1) 

ax 

where the constant a has the dimensions of volt cm and is given in terms 
of the electronic charge q and the dielectric constant k by 

4x0 
a = — - 

K 

For germanium a = 1.17 X 10"' volt em. 

The hole and electron flow densities Jp and /„ are" 



J„ = ^b{tiEn+D^]= ~hi.n 



E + 



g ax J 



where n and D = n kT/q are the hole mobility and diffusion constant re- 
spectively, k is Boltzmann's constant (8.63 X 10"^ ev per °C) and T is 
the absolute temperature. The ratio h of electron mobility to hole mo- 
bility we take to be unity. This makes the problem symmetrical in n and 
p and consequently easier to understand. Section VI will extend the re- 
sults to the general case of arbitrary h. 

Charge and Particle Flow 

For some purposes it helps to express the flow not in terms of Jp and 
J"„ but rather in terms of the current density I and the flow density 
y = Jp + J„ of particles, or carriers. The current density / = q{Jp — 
/„). Each carrier, hole or electron, gives a positive contribution to 
J if it goes in the +x direction and a negative contribution if it goes in 
the —X direction. In other words, / is the net flow of carriers regardless 
of their charge sign. The current / is constant throughout the intrinsic 

* See, for example, Electrons and Holes in Semiconductors, by W. Shockley. 
D. Van NoBtrand Co., New York, 1950. 
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region. Particle flow is away from the center of the intrinsic region. 
Carriers are generated in the intrinsic region and flow out at the two ends, 
the electrons going out on the A'' side and holes on the P side. Thus J is 
positive near the IP junction and negative near the NI junction. 
From the definitions of / and / and equations (2.2) 

- = ii-Eip + n) - D^ip- n) 
Q dx 

(2.3) 

J = nE(p -n) - D~ip-\-n) 
ax 

It is convenient to express the equations in terms of E and a dimen- 
sionless variable 

which measures how "swept" the region is. In normal intrinsic material 
s = 1. In a completely swept region s = 0; at the junctions with highly 
extrinsic material s ^ 1. Using Poisson's equation to express p ~ n in 
terms of E, equations (2.3) become 

/ = E ~ ^ — 

a dx^ 

where tr = 2 /i n,g is the conductivity of intrinsic material. The particle 
flow / is thus seen to be the gradient of a flow potential that depends 
only on E and s. 

Equations (2.5) can be written in the form 



(2.6) 

(2.7) 



r T, ■> d'E~\ 

where £ = \/kT/2an^ is the Debye length in intrinsic material and 
iJ.= 2y^=v|7: (2.8) 

is a field characteristic of the material and temperature. Specifically £, 
is \/2 times the field requii-ed to give a voltage drop kT/q in a Debye 



1246 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1956 

length. For germanium at room temperature £ = 6.8 10"^ cm and Ei = 
383 volts per cm. 

Both / and J are the sum of a drift term and a diffusion term. For 
charge neutrality, where p — nis small compared to p -\- n, both charge 
diffusion and particle drift can be neglected. We shall see later that, 
except right at the junctions, charge diffusion is negligible. 

The Equations of Continuity 
The two equations of continuity are 



dJ„ dJ„ 



V _ 



dx dx 



(2.9) 



where g is the rate of pair generation and r the rate of recombination. In 
terms of / and J, these become 



^ = (2.10) 

dx 



or / = constant and 



^ = 2(^ - r) (2.11) 

ax 

which says that the gradient of particle flow is equal to the net rate of 
particle generation, that is, twice the net rate of pair generation. 

To complete the statement of the problem it remains to express g and 
r in terms of n and p. 

Generation and Recombination 

The direct generation and recombination of holes and electrons follows 
the mass action law, in which g - r'ls proportional to ni - np. The con- 
stant of proportionality can be defined in terras of a Hfetime r as fol- 
lows: Let by = 5n « n; be a small disturbance in carrier density. Then 
defining t(j/ — r) = —Sn, we see that the proportionality constant in 
the mass action law is (2ntr)~\ So 



'='"'= 2n,-r 



nl - np ^2.12) 



and the generation rate 



2r 
is independent of carrier concentration. 



g = ^ (2.13) 
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In actual semiconducting materials, recombination is not direct. 
Rather it occurs through a trap, or recombination center. The statistics 
of indirect recombination has been treated by Shockley and Read^ for a 
recombination center having an arbitrary energy level Ei somewhere in 
the energy gap. At any temperature the trap le^'el can be expressed by 
the values ni and pi which n and p would have if, at that temperature, 
the Fermi level were at the trap level. Shockley and Read showed that, 
at a given temperature, the lifetime for small disturbances in carrier 
density is a maximum in intrinsic material. It drops to hmiting values 
T„o and Tpi, in highly extrinsic n and p materia], respectively. The foi-mula 
for {7 — r in terms of n and p is 

_ ^ n-i^ — np , , 

Tpo(n + ?li) + r^(p + Pi) 

For our purposes it is more convenient to define the lifetime t not by 
■^(ff — r) = —8n <C Ui , but rather as the proportionality factor in the 
mass action law. Then r is not necessarily constant independent of carrier 
density. From (2.12) and (2.14) 

Tpoin + ni) + Tnoip + Pi) ,^ , -^ 

' ^ 2^, ^2-1^^ 

We shall be interested in the lifetime in the region where n and p arc 
equal to or les-s tlian ??, . t decreases as n and p decrease; that is, t is less 
in a swept region than in normal intrinsic material. Let t = Ti for n = 
p = Ui and T = To for n = p = 0. The total range of variation of t is by 
a factor of 

'-^ = 1 + "--^"^ + "-"^ (2.16) 

TO PlTnO + UiTpO 

Let the energy levels be measured relative to the intrinsic level, and 
define a level So by 

In 4/^" 
Then if £( = So , riiTpo = PiT„(, . Now eq. (2.16) becomes 

Thus the variation in t increases as the ratio of t„o to Tpo deviates from 
unity and as the trap level moves away from the level fio ■ 

3 W. Shockley and W. T. Read, Jr., Phys. Rev., 87, p. 835, 1952. 



So = kT h 

TpO 
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The data of Burton, Hull, Morin, and Severien^ shows that a typical 
value of the ratio of r^o and t„o is about 10. Tliis means that the varia- 
tion ill T with carrier concentration will be less than 10 per cent provided 
&t is about 4/cT from So . In what follows we shall assume that this is so. 
Then we have the mass action law (2.12) with t a constant, which could 
be measured by one of the standard techniques involving small dis- 
turbances in carrier density. The general case of variable t is considered 
briefly at the end of Section W. 

Outline of the Solution 

To conclude this section, we discuss briefly the form of the equations 
and the solution in different parts of the mtrinsic region. First consider 
(2.6) for the current in the ideal case of equal mobilities. In Sections III 
and V we shall show that throughout almost all of the intrinsic region the ' 
current flows mainly by pure drift so we can take I = tsE. The reason 
for this is as follows. The quantity £^ is so small that the diffusion term 
remains negligible unless the second derivative of E becomes large — so 
large in fact that the E versus x curve bends sharply upward and both 
the drift and diffusion terms become large compared to the current I. 
This is the situation at the junction where / is the small difference be- 
tween large drift and diffusion terms. Thus (2.6) has two limitmg forms: 

(1) Except at the junctions the current is almost pure drift so / = 
asE is a good approximation. In Section III we derive an upper limit 
for the error introduced by this approximation and show how the 
approximate solution can be corrected to take account of the diffusion 

term. 

(2) At the junction, the drift term becomes important and the 
current rapidly becomes a small difference between its drift and diffusion 
terms and the solution approaches the zero current solution, for which 
sE = £" (fE/dx^. In Section V we derive an approximate solution that 
joins onto the / = usE solution near the junction and then turns con- 
tinously and rapidly into the zero current solution. We shall call this the 
junction solution. 

The abrupt change in the solution from (1) to (2) near the junction 
is shown to be related to a basic instability in the differential equation. 
This makes it impractical to solve the equations on a machine. 

When the applied bias is large compared to the built-in voltage drop, 
the junction region will be of relatively httle interest so the / = <tsE 
solution can be used throughout. 

In the region where 7 = csE there are two overlapping regions in 
which the equations assume a simple form. These are the following: 

* Burton, Hull, Morin and Severiens, J. Phys. Chem., 57, p. 853, 1953. 
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The No-Recombination Solution 

Here recombination is small compared to generation, r <^ g. This will 
be so in at least part of the intrinsic region for reverse biases of more 
than a few kT/g. The E versus x curve turns out to be given by a 
simple, cubic algebraic equation. 

The Recombination, or Charge Neutrality, Solution 

Here n — p is small compared to n + p, so the particle flow is by dif- 
fusion. We shall find that the 5 versus x curve is given by a third degree 
elliptic integral. As we move away from the center of the intrinsic 
region and toward the junctions, recombination becomes small com- 
pared to generation and the recombination solution goes into the no- 
recombination solution. In the region where both solutions hold, the 
solution has the simple form s = I/aE = A — x^ where A is a constant 
that must be less than | and the unit of length is twice the diffusion 
length. 

As the bias on an NIP structure is increased and the space charge 
penetrates through the intrinsic region, the region \\'here the recombina- 
tion is important will shrink and eventually disappear. 

Fig. 1 is a schematic plot of the field distributioii for the case where 
the applied bias is large compared to the built-in potential drop but not 
large enough to sweep all the carriers out of the intrinsic region. As the 
voltage is increased, the drop in field in the intrinsic region will become 
less and finally the field distribution will be almost flat from junction to 
junction. Only half of the intrinsic region is shown in Fig. 1. For equal 
mobilities the field distribution will be symmetrical about the center 
Xi of the intrinsic region. 

The illustration shows the recombination solution (1), which holds 
near the center of the intrinsic region and overlaps (2), the no-recom- 
bination solution. The junction solution (3) joins continuously onto the 
no-recombination solution at the point Xo and rapidly breaks away and 
approaches the zero-current solution at the junction. The figure is sche- 
matic and has not been drawn to scale. In most cases of interest, the low 
fields in the recombination region will be much lower and the junction 
solution will hold over a smaller fraction of the intrinsic region. 

It is convenient to take x = not at the center Xi of the intrinsic 
region but at the minimum on the no-recombination solution. As the 
applied bias increases, Xi approaches zero. 

Unegual Mobilities 

In the general case of unequal mobilities, it is no longer so that / is 
pure drift except at the jimctions. However we can define a linear com- 



1250 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1956 




Fig. 1 — Schematic of the field distribution and the three overlapping solutions. 

bination of Jp and /„ which has the same form as / in (2.6) and in 
which the diffusion term is negligible except near the junction. As we 
show in section VI, the effect of unequal mobilities is (1) to introduce 
some asymmetry into the curve in the region where the curvature is 
upward and (2) to displace the curve toward the NI junction (for the 
case where the electrons have the higher mobility). Thus the field is 
higher on the side w^here the carrier mobility is lower, as would be ex- 
pected. 

III. THE NO-RECOMBINATION CASE 

This section deals with the case where recombination can be neglected 
in comparison with generation. This will be so where n-p is small com- 
pared to Hi, 

The continuity equation for J now becomes 



dJ „ tu 
ax T 



(3.1) 
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Combining this with (2.7) gives 

rfi^ (^ " 7 " 2l? ^^'^^ 

where Li = Dt is the diffusion length in intrinsic material. 

Equation (3.2) can be immediately integrated. There are two con- 
stants of integration, one of which can be made to vanish by choosing 
3; = at the center of the intrinsic region, where the first derivatives of 
E and s vanish. {E is a minimum here and s a maximum). The solution 
obtained by two integrations is 

As we shall see later, the constant A is determined by the voltage drop 
across the unit. 

The exact procedure now would be to substitute 5 from (3.3) into (2.6). 
The resulting second order differential equation could, in principle, 
then be solved for E versus x. The exact solution, however, would be 
quite difficult. We shall discuss it in Section V. Here we make the assump- 
tion that throughout the intrinsic region the charge flow is mainly by 
drift, so that we can neglect the diffusion term in (2.6) and take / = 
a-sE, as discussed in Section IT. Later in this section we find an upper 
limit on the error due to this assumption and show how the cubic can be 
corrected to take account of the diffusion term. 

Putting s = I/<jE in (3.3) gives a cubic equation 



\eJ cE, 



for E/E\ as a function of x/2L{ . This equation contains two parameters 
/ and A. A determines the voltage and / is determined by the length 
2L of the intrinsic region. The relation is as follows: Let the applied 
voltage drop across each junction be at least a few JiT/g. Then the 
minority carrier currents from the extrinsic regions will have reached 
their saturation values. Call I^ the sum of the hole current from the 
A'^ region and the electron current from the P region. 7g comes from pairs 
generated in the extrinsic regions near the junctions. I, can be made 
arbitrarily small by making the N and P regions sufficiently highly 
doped (provided the diffusion length in the extrinsic material does not 
decrease with doping faster than the majority carrier concentration in- 
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creases). The current generated in the intrinsic region is qg per unit 
volume. So the density of current from pairs generated in the intrmsic 
layer is 2Lqg = qutL/T. Hence 

/ = /. + «^ 

T 

In what follows we shall assume that I, is negligibly small compared to 
7, Then 

Thus / is L/Li times a characteristic current equal to (1) the diffusion 
current produced by a gradient nt/Li or (2) the drift current produced by 
a field that gives the voltage drop kT/g in two diffusion lengths in normal 
intrinsic material. In germanium this characteristic current is about 5 
milliamperes per cm'. 

That the current I is proportional to L and independent of voltage 
follows from the neglect of recombination. When recombination is small 
compared to generation, then the current has reached its maximum, or 
saturation, value. All the carriers generated in the intrinsic region are 
swept out before recombining. It will sometimes be convenient to take 
<tEi as the unit of current. From the above and (2.8) 

_L = "V^^^ (3 5) 

In germanium <tEi is about 7 amperes per cm^ In general we will be deal- 
ing with currents that are small compared to this. For example, if L,- 
is 1 mm, we would have to sweep out an intrinsic region 3 meters long 
in order to get a current this large. If we take Ei as the unit field, cEi as 
the unit current and 2L. as the unit length then the cubic becomes E'^ - 
I/E ^ x^ - A. 

For a given structure and temperature the field versus x curves form 
a one parameter family. A determines both the field distribution and the 
voltage. The voltage increases as A decreases. Fig. 2 is a plot of E/Ei 
versus x/2Li for L/2Li = 0.1 and several different values of A. Fig. 3 
is for L = 2Li and Fig. 4 for L/2Li = 3. 

There is an upper limit on A but not lower limit. The reason is as 
follows: As A increases, the minimum value of E (at x = 0) decreases 
and the maxunum value of s increases. So if A is too large, the maxi- 
mum s will be so large that we cannot neglect recombination, %vhich 
becomes important when np approaches n^^ or s approaches 1. Fre- 
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quently recombination can bo neglected over parts of the intrinsic 
region but not near the center, where the field is a minimum and the car- 
rier concentration a maximum. Then (3.4) will represent the field chs- 
tribution over that part of the region where recombination is unimpor- 
tant. The correct solution will join onto the cubic as we move away 
from the center of the intrinsic region, which will no longer be at the 
X = point on the cubic. In Section IV ^e solve the equations for the 
recombination region and show how the solution approaches the cubic. 
We will show that, as A increases, the distance from the center of the 
intrinsic region to the x = point on the cubic also increases. The 
value A = ^ corresponds to an infinitely long intrinsic region. For a 
larger A there exists no exact solution that could join onto the cubic 
as recombination becomes negligible. In Figs. 3 and 4 the A = ^ curves 
join onto recombination solutions at values of E which are too low to 
show. 
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Fig. 2 — Field Distributions for L = 0.2Li . 
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As A decreases and becomes negative the cubic approaches the form 

^i ^ \ (3.6) 



E' = E,' + £i^ 



2Li 



where E<i^ = ~AEi is the minimum value of £^^ This form of the solu- 
tion will be valid when the minimum E is large compared to {lEi /a). 
As Eo increases, the voltage increases and the curve becomes flatter. 
This is because the increasing field sweeps the carriers out and reduces 
the space charge; so the drop in field decreases. 

If (3.4) for E/Ei versus .r/2L.- is extended to indefinitely large values 
of x/2Li , it approaches the straight line of slope 1 going through the 
origin. Since E is always positive the curve is above this straight line at 
X = 0. If ^ is negative the cui-ve is always above the straight line and 
always concave upward. If A is positive, the curve crosses the straight 



1.5 








































y 


/ 
















y 


y 












A = -) 


^ 


^ 


^ 














^ 


■^ 


































^ 




















// 


















A 


/ 


/ 














/' 


/ 


/ 
















/ 


/o.i / 


/ 


/ 












/ 


v 




/■/3 


/ 










/ 


/ 




/ 




1 










/ 




/ 




/ 




^A=2/3 




. 


^ 




/ 




J 






/ 









___ 


^ 




. 


J 




^ 


' 





01 2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
X/2Lt 



Fig. 3 — Field Distributions for h = 2L( . 
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line at B/Ei = I/cEiA and thereafter remains under it approaching it 
from below. For positive A the curvature, which is upward near the ori- 
gin, changes to downward at about x/2Li = \fA. 

The carrier concentrations n and -p can be found from the E versus x 
curves with the aid of Poisson's equation p — 71 = 1/a dE/dx and the 
definition s = (n + p)/2n,- with s = I/aE. These relations and (3.4) 
give 



(3.7) 




^ 2aEy 



From (3.4) and (3.7) we may distinguish the following two regions 
on the cubic: 

(1) When E^/Ei is smaller than I/aEi (which as we have seen is 
usually smaller than unity), the E versus x curve is concave upward, the 
hole and electron concentrations are almost equal (charge neutraHty) 
and the particle flow is by diffusion. 

(2) When E^/E^ is greater than I/<tEi , in general there is space 
charge and the particle flow, like the charge flow, is by drift. The curve 
is concave downward for positive A. 

Figure 6, which we will discuss in Section IV, shows the field and car- 
rier distributions for L = 2L, and A = 0.665 plotted on a logarithmic 
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scale to show the behavior at low values of field and carrier density. 
In the region of no-recombination the field distribution is indistinguish- 
able from that for A = I, which is plotted in Fig. 3 on a Imear scale. In 
the region where recombination is important the sohition is found from 
the assumption of charge neutraUty as will be discussed m Section IV. 
The cubic and charge neutrality solutions are each shown dashed outside 
of their respective ranges of validity. For A = 0.665 the half length of 
the intrinsic region is 2.098 X 2L, . Thus the length of the intrinsic re- 
gion is more than twice the effective length 2L m which current is 
generated. The effective length will be discussed in more detail in Section 
IV and it will be shown that the effective length 2L of current generation 
is equal to the twice the distance from the IP junction to the minimum on 
the cubic. As explained earlier, it is convenient to take a; = at the mini- 
mum on the cubic. 

Inirinsic-Exirinsic Junction Under Large Bias 

Consider the limiting case of an intrinsic-extrinsic junction as the 
bias is increased and the space charge penetrates many diffusion lengths 
into the intrinsic material. Then the field distribution approaches the 
straight fine E/Ei = x/2Li . This, by Poisson's equation, means that 
there is a constant charge density of Ni where 

AT = ^ = V2£^ 
2aLi Li 

Thus in the hmit, the field in the intrinsic region approaches that in a 
completely swept extrinsic i-cgion havmg a fixed charge density of A^ . 
In germanium at room temperature Ni is about 4.10'" cm"^ As the 
field approaches the limiting form, the voltage V approaches EiL /4L, . 
Thus the limiting form of the current voltage curve is 



cE 



\ L. r 2E,L 



So in the limit the current varies as the square root of the voltage. Typical 
values for germanium at room temperature are aEi = 7 amps cm , 
£/Li = 10"^ and 2EiLi = 50 volts. 

Equivalent Generation Length for an Intrinsic-Extrinsic Junction 

It should be noted that for an IP structure the current is the same as 
for an NIP structm-e with an infinite 7 region, or at least an I region 
that is long compared to the distance of penetration of the space charge. 
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Thus the equivalent length of current generation is 2L even thoutrh the 
current is actually being generated in an effective length L. The reason 
is that for an NIP structure the holes entering the right hand half of the 
I region were generated in the left hand side. For an IP structure the 
holes entering the space charge regions from the left were injected at the 
external left hand contact to the I region. 

Applied Voltage 

In all cases the voltage can be found from the area under the E versus 
X curve. In Figs. 2 to 4 the area under the curves gives the voltage ac- 
curately; recombination becomes important only where the field is so 
low as to have a negligible effect on the total voltage drop. 

Correction oj the Cubic 

To conclude this section we consider the error introduced by using the 
assumption I = asE. For simplicity take Ei as the unit field, 2Li as 
the unit length and aEi as the unit current. Then the cubic becomes 
E^ — I/E = x^ ~ A. The corresponding exact solution is E^ ~ s = 
x^ — A where the relation between s and E is given by equation (2.6) 
which in dimensionless form is 

.'§ = .!.-/ (3.8) 

where £^ is of the order of 10~^. 

Let 5E and Ss represent the difference between the cubic and the cor- 
rect solution at a given x. Assume that SE and its second derivative are 
small compared to E and its second derivative respectively. Then 5s = 
2E8E and on the correct solution sE - I = (2^' + I/E)BE. So (3.8) 
becomes 

E \2E' + l) dx' "-"^'^^ 

To obtain a first approximation to 5E/E we use the cubic to evaluate 
^Ejdx. It is convenient to express the results in terms of a dimension- 
loss variable z = E/I^'^, or if E and / are measured in conventional units 

z = E{a/Exl)"\ Then (3.9) becomes 

hE 1/L,£VV ^ Y , (x\z\\~z') 



+ [^]^^r^. (3.10) 



E 2 V LV \3' + 1/ \2L/ {z^ + I)* 
fvhen the lengths are in conventional units. 
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The first term has a maximum value of 0.35 (L,£/L^)^'^ at z = 0.6 
and the second term a maximmn value of 0.18 atz = 0.5 and x = L. 

The dashed curve in Fig. 2 for A = .01 is the corrected cubic. For 
the other curves in Fig. 2, the correction is smaller. For the curves in 
Figs. 3 and 4 the correction is too small to show. 

Limits on the Solution 

We now show that SE as derived above is not only a first approxima- 
tion but also upper limit on the correction necessary to take account of 
charge diffusion. That is, an exact solution to (3.8) lies between the cubic 
and the corrected cubic. 

Consider the region where the second derivative of E is positive so 
that the perturbed curve lies above the cubic as in Fig. 2. On the cubic 
we have sE ~ I = 0. As we move upward from the cubic and toward the 
dashed curve, sE - I increases. The value oi sE ~ I on the dashed curve 
just equals the value of Z^ cfE/dx^ on the cubic. However, the dashed 
curve has a smaller second derivative than the cubic. Thus in moving 
upward from the cubic toward the dashed curve sE ~ I increases from 
zero and £.^(fE/dx^, which is positive, decreases; on the dashed curve 
sE — I is actually greater. Therefore there is a curve lying just under 
the dashed cui-ve where the two sides of (3.8) are equal. The same argu- 
ment applied to the region where the second derivative is negative shows 
that the equation is satisfied by a curve lying just above the first per- 
turbation of the cubic. Where the curvature changes sign, the cubic is 
correct. 

It should be emphasized again that the neglect of the diffusion term 
in the current is justified only for the ideal case of equal hole and electron 
mobilities. For unequal mobilities both drift and diffusion will be im- 
portant in current fiow. However, as we will discuss in section 5, we can 
simplify the problem of unequal mobilities by defining a fictitious current 
that has the same form as I in (2.6) and (3.8). 

IV. RECOMBINATION 

As discussed in Section III, when the voltage for a given current is re- 
duced, s increases and near .r = becomes comparable to unity. Then 
recombination becomes important and the cubic solution breaks down, 
or rather joms onto a solution that takes account of recombination. 
When recombination is important the center Xi of the intrinsic region is 
no longer at the x = point on the cubic but to the left of it. That is, if 
we want the same current with continually decreasing voltage, we even- 
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tually get to the point where a longer intrinsic region is required. Finally 
for a given current we reach a minimum voltage which corresponds to an 
infinite length of intrinsic region. Another way of saying this is that, 
when recombination becomes important, the length L defined in terms of 
the current by / = qg2L = qni/rL is no longer the half length of the 
intrinsic region. 

Equivalent Generation Length 

We shall continue to define Z, by / = qriifrL. Thus h is an equivalent, 
or effective, half length of current generation and not the half length of 
the intrinsic region. By definition L is the length such that the amount 
of generation alone in the length L is equal to the net amount of genera- 
tion (generation minus recombination) in the total half length of the 
intrinsic region. Hence 

qL= f " (g - r)dx (4.1) 

where a;, is at the center of the intrinsic region and Xp at the IP jimction. 
We shall for the most part deal with reverse biases of at least a few kT/g, 
in which case recombination is negligible at the junctions. Then the exact 
solution becomes the no-recombination solution before reaching the junc- 
tions. Wc shall continue to take x = at the puiiit ilE/dx = ds/dx = 
on the no-recombination solution which the exact solution approaches 
as recombination becomes negligible. 

Simplifying Assumptions 

The general differential equation with recombination will be the 
same as for no-recombination except that g — r replaces g. From (3.1) 
and (3.2) 

£(l-)=,i{-:) 

From (2.12) and (2.13) and Poisson's equation 

r _np _ fn + pY (n - p)^ ^ ^, ^ ^ f^^Y (43) 



g iif \ 2ni / (2?i,) V-^i dx 

The following analysis will be based on the assumption of charge neu- 
trality. That is we neglect terms in p — n in comparison with those in 
p + 71. In particular this means: 

(1) The charge flows by drift so / ^ asE. This is the same assumption 
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made in the no-recombination case. It will be an even better approxima- 
tion in the recombination region, where the second derivative of E is less. 

(2) The particle flow is by diffusion. That is, E^/Er can be neglected 
in comparison with s. 

(3) The ratio of recombination rate r to generation rate g is propor- 
tional tog — r\ that is gr — r = ff(l — s ). 

All of these simplifying assumptions can be justified by substituting 
the resulting solution into the original expressions and showing that the 
neglected terms are small when recombination is important. If the 
solution is substituted into (4.3) and (2.6) the neglected terms wall 
turn out to be negligible — and therefore assumptions (1) and (3), 
justified — when s^ is large compared to £/L. . Assumption (2) follows 
from (1) and the fact that I/aEi is small compared to unity. 

Assumptions (2) and (3) may also be justified by the discussion fol- 
lowing (3.7) in the following way: Where recombination is important s 
must be near unity. So the cubic will begin to break down when s = 
I/<jE becomes near to unity, or when E approaches I/c. However, if 
E is approximately I /a then uE^jJEx is approxunately {I/aEif, which, 
as we saw in the Section III, is small compared to unity in practical 
cases. Thus recombination becomes important and the solution joins 
onto the cubic in the range w^here E^/Ei is small compared to I/aEi . 
In this range the particle flow is by diffusion and p - n is small compared 
to p -h n. As wc move toward the center of the intrinsic region s increases 
and E and dE/dx decrease. Therefore, since assumptions (2) and (3) 
are good where the solution joins onto the cubic, they are good through- 
out the region where recombination is important. 



The Recombination Solution 
The diiTerential equation (4.2) now takes the form 

(4.4) 



dh (1 - s') 



The solution for s in the recombination range is seen to be the same for 
all values of the current. When s has been found E is found from E = 

I/<JS. 

For small disturbances in normal carrier concentration, s is only 
slightly different from unity and (4.4) takes the famiflar form 
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which says that the disturbance in carrier concentration varies expo- 
nentially as xjhi . 

Equation (4.4) can be integrated once to give 

where So is the value of s at the center of the intrinsic region where s is a 
maximum. 

As recombination becomes unimportant, s^ becomes small compared 
to unity and (4.5) approaches the form 

and the solution joins onto the no- recombination solution. 

Joining onto the Cvbic. 

We have seen that the solution joins onto the no recombination solu- 
tion, in the region where particle flow is by diffusion so that the no 
recombination solution has the form s = A ~ (.T/2L,■)^ This is readily 
transformed to the form (4.6) with 

A = so (l - ^) (4.7) 

Thus the one arbitrary parameter So in the recombination solution 
determines the parameter A in the cubic that the recombination solution 
approaches. Since the maximum value of So mider reverse bias is unity, 
the maximum value of A is f . Negative values of .4 correspond to solu- 
tions where recombination is always neghgible. 

The s versus x Curve 

To find s versus x we integrate (4.5). There is one constant of integra- 
tion, which is fbced by the choice of a; = 0. We have taken a; = at the 
point where dE/dx = ds/dx = on the cubic. To make the recombina- 
tion solution join the cubic we choose the constant of integration so that 
the recombination solution extrapolates to s = at (a:/2L,)^ = A . Then 

_L = /a -^^ [' t^s 

2Li ^ 2 io V3(so - s) " (so3 - s3) ^ -^^ 

which can be expressed in terms of elliptic integrals. 
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Fig. 5 — Variation of s = p/nt = n/tii in the range where recombination is 
important. 

Deep in an infinitely long intrinsic region the carrier densities ap- 
proach their normal values n = p = n, , or s = 1. Putting So = 1 in 
(4.8), we find that as s approaches so ^ 1, x becomes infinite. This will 
be the solution for a simple intrinsic-extrinsic junction. Fig. 5 is a plot 
of s versus x for various values of s,, . Tlie dashed curves represent the 
corresponding no-recombination solution s = A ~ (.r/2Lv) . 

The IP Junction 

It remains to find the position of the IP boundary. We now show that 
if recombination is unimportant at the junction, so that the solution 
joins onto a no-recombination solution, then the position of the junction 
is at X = L where L is the effective length of current generation and 
X = is the point where dE/dx = ds/dx = on the no-recombination 
solution (which of course will not be valid at x = 0). The proof is as 
follows: From the definition (4.1) of L and (4.2) 



d_ 
dx 



— s 



(4.9) 



r*p >. f''^ d^ / E^ 

j^ = I (1 - r/,) d. = 2Lf I I, (^ - 

= 2L,^ [ 

If the boundary comes where recombination is negligible so that 
{E/Eif - s = (x/2L,)' - A, then (4.9) gives Xp - L. Physically 



dx 
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this means that the amount of recombination in the interval from 
a; = to a; = L is just equal to the excess amount of generation in the 
interval from the center of the intrinsic region to x = 0. 

If the applied reverse I)ias is less than a few kT/rj then recombination is 
important even at the junction and there is no joining onto a no-recom- 
bination solution. In this case (4.9) says that for a given choice of cur- 
rent (and hence of L) the boundary comes where 

Example. Fig. 6, which we discussed briefly in Section III, is a plot 
of the field and carrier distributions for L = 2Li and So = 0.95, for which 
A = 0.665. The hole and electron densities were found from (3.7) and 
p + n = 2n.s where s is found from Fig. 5. When s approaches So (4.8) 
for X versus s takes the simple form 

This will be accurate when Sn — s is small compared to 1/so — So . We 
have used (4.11) to evaluate the s versus x curve beyond the range of 
the So = 0.95 curve in Fig. 5. 

It is seen that the recombination solution in Fig. 6 joins the cubic 
in the range where n and p arc still almost equal. 

Variable Lifetime, 

Finally consider the general case where the variation in t with car- 
rier density cannot be ncglectetl. Then, with n = p = nis, (2.15) be- 
comes t = To + (t,- — T[,)sand//,'in (4.4) is replaced by /)r[l + (t,/to — 
l)s] where n/To is given by (2.17). The more general form of (4.4) can 
be solved graphically after one integration. The solution will join onto a 
cubic if (ti/to — l)s becomes small compared to unity before space 
charge becomes important. This will be so if (t,/to — if^I/aEi is small 
compared to unity. 

V. THE JUNCTION SOLUTION 

In this section we consider the solution near the junctions, where 
the assiunption / = csE breaks down. We shall deal with reverse biases 
of at least a few hT/q so that recombination is negligible at the junctions. 

The junction solution will therefore join onto tlie no-recombination 

solution. We shall use the cubic solution in the no recombination region. 

Again it is convenient to use dimensionless variables with Ei as the 
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Fig. 6 — Field and carrier distributiona for L = 2i,i and A = 0.665 (so = 0.95). 

unit field, 2Li as unit length and aEi as unit current. Then on the cubic 
s = I/E, and E^ — I/E = x'^ ~ A. The current is related to L by 
7 = ■\/2£L where the dimensionless £ is of the order of 10"' for ger- 
manium at room temperature. Substituting the exact no-recombination 
solution E- - s = x^ ~ A into the solution (2.6), or (3.8), for the current 



gives the second order differential equation 



<fE 



(5.1) 



for £ as a function of x. The two boundary conditions are as follows: At 
X = 0, dE/dx = by symmetry. At the IP junction the carrier concen- 
tration must rise and approach that in the normal P material. For u 
strongly extrinsic P region the normal hole concentration P is large 
compared to both n, and the electron concentration. Thus s must in- 
crease and approach P/2;?, » 1 as we approach the P region. Clearly 
the cubic cannot satisfy this requirement. On the cubic the maximum 
value of s comes at x = and is less than unity. As we approach the June- 
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tion E increases so 5 = I/E must decrease. Thus the correct solution 
must break away from the cubic near the junction. 

Instability of the Solution 

Equation (5.1) has two limiting forms and makes a rather abrupt 
transition bet^veen them. Over mo.st of the intrinsic region, the quantity 
in brackets [Es — /] = [E{E^ — a:" + A) — I] almost vanishes. It differs 
from zero just enough that when multiplied by the very large factor 
£^ fn 10" it gives the correct second deri\•ati^'e of E. In Section III we 
derived an upper limit on the small deviation 5E from the cubic required 
to satisfy the differential equation. If E deviates from the cubic by more 
than this small amount, then the second derivative of E becomes too 
large. This increases the deviation from the cubic, which further in- 
creases the second derivative and so on. E and s rapidly approach 
infinity in a short distance. This, of course, is the required behavior at 
the junction. The rapid increase in s makes it possible for s to approach 
P/2ni . 

In Section III we showed that there is a solution to the differential 
equation that lies within a small interval 5E from the cubic. Suppose we 
try to solve (5.1) graphically or on a machine starting at x = 0. There 
are two boundary conditions: By symmetry dE/dx = at 3; = 0. We 
choose for £"(0) a value somewhere in the interval dE{Q). The resulting 
solution will not long remain in the interval 5E(x). In fact there is only 
one choice of E(0) for which the solution remains close to the cubic 
from .X- = to X = cc . For any other E{0) the solution would remain 
close to the cubic for a certain distance and then abruptly become un- 
stable and both E and s approach infinity. £"(0) must be so chosen that 
the solution becomes imstable and E and s become large at the junction. 
However it is impractical to set E(0) on a machine with sufficient ac- 
curacy to insure that the solution will remain stable for a reasonable 
distance. A more practical procedure is to find a solution which holds 
near the junction and joins the cubic to a solution in the adjacent 
extrinsic region. 

Zero Bias 

It may be helpful to approach the junction solution by reviewing the 
simple case of an IP junction under zero bias. Both charge and particle 
flow vanish. The vanishing of particle flow means that in the intrinsic 
region E^ — s is constant, (2.7). Since E = and s = 1 in the normal 
intrinsic material, it follows that E' — s = 1. With / = the equation 
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for current becomes 

This can be integrated at once. The boundary conditions are dE/dx = 
when E ^ and E = Ej at x = L; the field E,- at the IP junction will 
be determined by joining the solutions for the 7 and P regions. The 
solution can best be expressed by parametric equations giving x and the 
potential V as functions of E. 

c^i dE r ._, E ,Ei 1 

y^-y = ^ i vi + E^/2 = T L""' VI " "^^ vd ^'-'^ 

where we have used the relation between dimensionless quantities 
£ = \/2kT/q, which follows from (2.8) with Ei = 1. It will be more 
convenient to express voltages in terms of kT/q rather than in terms of 
the unit voltage 2EiL, ; then the ratio qV/kT is independent of the 
units. For convenience we take the voltage as increasing in going toward 
the IP jmiction with 7 = in the normal P material. The voltage Fj 
at the junction is found by joining solutions. 

On the P side, let the excess acceptor density be P. Adding the term 
— aP to the right hand side of Poisson's (2.1), and proceeding as before 
we have, instead of (2.5) 

where Sp - P/2ni . We shall assume that the P region is strongly ex- 
trinsic so that n « p. Then s = Sp in the normal p material, where 
E = V = 0. Hence 

;.'-. = .,(g-i) (5.5) 

In the intrinsic material the corresponding solution is E^ - s = -1. 
Since both E and s are continuous at the junction, qVj/kT = I — 1/sp 
where 1/Sj, can be neglected. Thus S/ = Sy = Sp exp [~(qVj/kT)] = 
Sp/e where e = 2.72 is the base of the natural logarithms. 

Knowing Ej we can find the field and potential distributions in the 
intrinsic material from (5.3) and (5.4). 
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Reverse Bias 

Now in the intriusic region, E' — s = x' — A. Let Ec be the value 
of E at the junction as given by the cubic, and let Sc = !/Ec be the 
corresponding value of s. Then at the junction x"^ — A = E^ — s^ . In 
the P material equation (5.5) will still be a good approximation near 
the junction, where the additional terms arising from the flow will be 
negligible. Joining the solutions for the / and P regions and neglecting 
Sc in comparison with Sp gives 

^-1^ = 1+^ 
hT ^ Sp 

Again using Sj = Sp exp [— CjVj/fcT')] we have 

Ef = E: + Sp exp [-(1 + £,V5p)] (5.6) 

In most practical cases £/ will be small compared to Sp = P/2n,- so 
Ej will be the same as for zero bias. 

Junction Solution 

We now consider an approximate solution that joins smoothly onto 
the cubic and has the required behavior at the junction. Let x = xq 
be the point where this solution is to join the cubic. Then in (5.1) x^ 
must lie between .Tq and L . We can obtain two limiting forms of the 
solution by giving x the two constant values, Xo and L respectively. 
It will be best to take x = .To since in practical cases the x^ term is not 
important except near the point where the junction solution joins the 
cubic. In all cases the uncertainty due to taking x~ = constant can be 
estimated by comparing the solutions for x = .ro and x = L. 

With X constant, (5.1) can easily be integrated. The two boundary 
conditions are (a) E = Ej at x = L, where Ej is given by (5.6), and (b) 
to insure a smooth joining, the slope at a; = .To must be the same as that 
of the cubic, namely 

/dE\ ^ 2X0 ,„. 

\dxJo 2Eo + I/Eo' ' ^ ' 



The first integration of (5.1) with .t = To gives 

where (dE/dx) is given by (5.7) and Eq^ — I/Eo = Xa — A. The E versus 



(5.8) 
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X curve can now be found from (5.8) and 

•'- ^''^ ^ (5.9) 

In general we will be intrested in cases where the junction solution 
holds over a length L — Xo that is small compared to L, so ^^-e can take 
Xo = L in (5.7). It will also be valid to let Eo in (5.7) and (5.8) be the 
value Ec on the cubic at x = L. Putting Ec = Ea in equation (5.6) then 
gives Ej in terms of E^ and Sp = P/27i,- , where P is the majority carrier 
concentration in the extrinsic region. In what follows we shah use these 
approximations. It will be convenient to express xa = L in (5.7) in 
terms of / using / = s/^L^. We continue to use dimensionless quan- 
tities with Ex , ILi and aEi as the units of field, length and current re- 
spectively, and 2L,^i as the unit of voltage. In general however we can 
express voltages in terms of kTjq. 

When E^ is either large or small compared to /, the junction solution 
takes a simple form and the field and potential distributions can be 
found analytically. We next consider two approximations that hold in 
those two cases respectively. Relatively good agreement between the 
two solutions at E^ = I indicates that each solution may be used up to 
this point. 



-1/2 

dE (5.10) 



Case of Ea Large Compared to I 
From (5.7) to (5.9) 

This can be solved in the following two overlapping ranges where the 
integrand has a simple form: 
Range 1. Here E - Eo'is small compared to 2Eo , so (5.10) becomes 



X — Xo = — 7^^ sinh 



I'D 



(K-E3)?|i] (5.11) 



V2Ea 

Since E and Eo are ahnost equal, we have for the voltage drop in this 
range 

y - Fo = Eoix - Xo) (5.12) 

Range 2. Here E^ - Eo is large compared to 2(£.L/Eof, so (5.10) 
gives 
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- - ^ = V2. /;,^ = ^ (ctn.- I - etnh- |) (5.13) 

From £'0^ » / it follows that Ranges 1 and 2 overlap. By joining the 
two solutions in the overlap region, the solution in Range 2 can be 
written as 

Putting £^ = Ej in (5.14) gives the distance over which the junction 
solution holds. In general we will be interested in cases where Ej is 
large compared to Eo so (5.14) becomes 

L ~ Xo 3tn2za ,. ,^. 

—[- = 2 -IT ^^-^^^ 

where I = s/^^/I^'^ and as before 2o = Ea/P'^. In conventional units 

Z = 2L (^^J" (5.16) 

Fig. 7 is a plot of {L — xo)/l versus zo . In germanium at room tempera- 
ture £Li will be around 10"^ cm. Thus the junction solution will hold 
over a region that is small compared to L if L is large compared to 
3 X 10~^ cm. 

Again it is convenient to use the relation £ = \/2kT/q to express 
the voltage in terms of kT/g. 

By joining the two solutions in the overlap region, the voltage in Range 
2 can be expressed as 

V-Vo = ^en (^-^^ (E' - LV) (5.18) 

Setting V = Vj and E = Ej in (5.18) gives the total voltage drop in 
the region where the juuction solution holds. Let AV bo the difference 
between Vj — Vo and the built in voltage drop at the junction. Then 
substituting (5.6) with Ec = Eg into (5.18) and subtracting the built in 
drop we have for AV, 



,. IcT . Eo Eo' /r in\ 

Al = — hi— (5.19) 
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I/Eo is equal to the I'aluc of s on the cubic at x = L. For positive values 
of A the maximum value of E^/I is L/I = l/v^£ as can be seen from 
the cubic. In germanium at room temperature £ is about 10 (for 
2Li = unit length) so the reverse bias produces an additional voltage 
drop in the junction region equal to about 7kT/g. For negative values 
of A the additional voltage drop near the junction would be higher. 

Comparing (5.3) and (5.13) we see that the junction solution reduces 
to the zero bias solution when E'' is large compared to Eo^ + 2. In this 
case both solutions have the simple form 



and 



L - 0. = V2£ (i - ^ 



.,.-.^?f.4 



(5.20) 



(5.21) 



Case of Eo Small Compared to I 

Now from (5.7) and (5.8) with xq ^ L = 7\/2£ we have 



dx/o \ I 



«■(?-!)] 



(5.22) 
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Again there are two overljipping ranges where the solution has a simple 
form : 

Range 1. Here E' is small compared to 2//£'o . This will be so even 
when E becomes large compared to Eu . Setting o" ^ '^.Ei/l aiid y = 
E ~- i'o in equation {5.22) and integrating gives 



X - x„ - £A/^ 



En i^~^' dt) 



= ii ^/^sinh 



fin . ■ I -1 / -^ ~ -tio 



C\ 



(o.23) 



and : 

V - 7o 

= Y \/^f (Vc,-^ + (E - E.r- - fi) + 2/i'o(.r - .r„) 

Range 2. Here E is large compared to Ei, . It follows from Ea « / 
that E is also large compai-ed to Ci . Setting c-r = 21 /En we have 

clE 
E ^/W^c> 

^ ( esch"' - - cs<'h- ^ 



L - X = V2£ / 



-l/f(^ 



c-i c_ 

Joining (5.21) and (5.23) whci'c they o\'ei'lap we have in I'ange (2 

E 



a: - :r„ = £ 4/±^fn{-^2 



r '" V^.V U + V^-T J ^-'-^^^ 



Putting .1- = L and /:' = A'; in (5.2r)) gives the lengtli L — Xt, in which 
the junction sokition holds. If Ej is large compared to cj , then 

^-' — ^<i , Ao ,4 /- „_v 

where as before 2o = E,,/I^''^ and ^ is given by (5.10). Fig. 7 is a plot of 
(L — Xo)/l versus za . The two approximations (5.15) and (5.27) for 
z„ « 1 and 2(1 » I respectively are shown dashed. Both become inaccu- 
rate as they are extended toward zt, ^ 1. The point at zo = 1 was ob- 
tained graphically. Each approximation is in error by about 28 per cent 
h(!re. The error will decrease as each approximation is extended away 
from Zi) = 1 toward its range of validity. 
The voltage in Range 2 is given by 



1272 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1956 

V, _ 7 = ?^ fsinh-' ^ - sinh-^ -] (5.28) 

Q L C2 C.J 

or again joining (5.28) to the solutions in Range 1 we liave in Range 2 

V -Vo = ^ sinh"^ - + 2Eo(x - xo) (5.29) 

q C2 

The total voltage drop in the junction can be found by setting V = Vj 
and E - Ej in (5.29). The term 2Eo(L - .ro) will be negligible. When 
£" is large compared to d + 2 the junction solution reduces to the 
zero current solution as can be seen by comparing (5.3) and (5.25). 
Then the solution has the simple form (5.20) and (5.21). Ej will always 
be large compared to Cg . (-E/ is approximately sje and ci = 2.So where 
So is the value of s where the junction solution joins the cubic.) Thus the 
difference AF between V, -~ Vn and the built in voltage is 

^V = --ln^ (5.30) 

Example. Fig. 8 shows the field distribution near the IP junction 
for the case L - 2Li and A = |, for which the intrinsic region is in- 
finitely long. The field distribution near the junction, however, will be 
indistinguishable from that for A = 0.665, or So = 0.95, for which the 
mtrinsic region is about twice the effective length of current generation. 
We have taken Ej ^ 30, which corresponds to an excess acceptor den- 
sity P = 4.7 X 10^ Ui in the P region. Over the range where the junc- 
tion solution holds the cubic gives an almost constant field E = En — E^. . 
The junction sokition goes from the cubic to the zero bias solution in a 
distance of the order of the Dcbye length. The sum of the built in volt- 
age and the voltage derived from the cubic differ from the correct voltage 
by the order of &Ei or about kT/q. 1"he total voltage is about 0.3 EiLi , 
which would be about 11 volts in germanium at room temperature. 

VI. GENERAL CASE, UNEQUAL MOBILITIES 

This Section deals with the general case where the ratio of the hole 
and electron mobilities is arbitrary. The procedure is similar to that 
used in the preceding Sections. Many of the results for & = 1 are useful 
in the present, more general, case. We shall deal first with the no-reeom- 
bination ease and again find that E is given by a cubic. However, the 
field distribution is no longer symmetrical and the coefficient of the I/E 
term in the cubic is a linear function of x instead of a constant. The 
differential equation for s in the recombination region remains un- 
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Fig. 8 — Field Distribution near the IP Junction for L = 2Li and A = 5. 

changed. It is no longer so that charge diffusion can be neglected except 
near the junctions. However, there is a Unear combination of Jp and J„ 
in which the diffusion term is negligible except near the junctions. 



Basic Relations 

The equations arc the two continuity (2.9) and Poisson's (2.1). The 
formulas for (/ — r remain unchanged, since they involve only the 
statistics of recombination and are independent of inobiHty. The hole 
and electron currents are given by (2.2) with b arbitrary. Equation (2.2) 
for Jp in terms of E, p and ?i remains unchanged. Now J„/b has the same 
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i'nnn ;is ./„ liad Tor the /) = 1 vane. It is therefore desirahle tu deal with 
(he fifiitioiiH carrier How Jp -\- J„/b and the fictitious current q(J,, — 
J ,Jb) since these have the same form in terms of E and s = {n -\- p)/2n; 
as ./ and / had for 6=1. Thus 



f/l -^p - -r ) = 



bj 1 + ?^ " L d^. 



(6.2) 



where Ei and £ lia^'e the same meaning as before and the conductivity 
of intrinsic material is now a ^ 9"im(1 + f')- -A-s before D and m are re- 
spectively the diffusion constant and mobility for holes. Equations (6.1) 
and (6.2) i-educe re.specti\-ely to (2.7) for ./ and (2.(i) for / = g(Jj. - J,.) 
where b = I. 

When the flow is by pure diffusion, the holes and electrons diffuse "in 
parallel" so (he effecti\'e (Mffusion constant is the reciprocal of the average 
of the ro(NprocaI hole and electron diffusion constants. Hence the effective 
diffusion [ength i.s given by 

Lf = Dr-^ (().3) 

1 -\- b 

We continue to let 2L = I/gg be the effective length of current genera 
tion; again it is the actual length for the no re(!oinbination case. Let x„ 
and Xp be the coordinates of the NI and //■* junctions respectively. 
Since the problem is not symmetrical we will not take j; = in the center 
of the intrinsic^ region even for the no- recombination case. 

No-Recombinalion Case 

Setting r = we can immediately integrate the continuity equations 

(Up ^ 
d.v 

subject to the Ijoundary conditions: 

at the NI junction, :r = x„ , 

(0.4) 
at the /P junction, ;f = Xp , 

The result is Jp ^ g{x — x„) and J„ ^ g{x — Xp). This agrees with / = 
q('^p ~ -fn) = 29(/L since 2L = Xp — x„ is the length of the intrinsic 
region^ which, for no-recombination, is also the effective length of cur- 



dJ„ 
dx - ' 




J„=o, 


J. = ~I/q 


Jv - I/q, 


Jn= 
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rent generation. It will be convenient to choose a: = so that x„ = 
~Xp/b. Then the origin is nearer to the NI junction for h > 1. Now 
from this and the boundary conditions (6.4) and / = 2ggL we ha^'^e the 
positions of the junctions: 

x„ 26 x„ -2 ^gg^ 



L 1 + 6' L 1 + 6 

As before, the junctions are at .-c = ± -6 for 6 ^ 1. 

We can now find the fictitious carrier flow Jp + J„/b and the fictitious 
current q{Jp — Jn/b) as functions of x. 

j^ + 4. = (l^) ,. (6.6) 



6 V 6 



('--*) 



^' ^+^^1 (6.7) 



1 + 6 V L 



where the dimensionless parameter )3 = (6 — l)/46. Thus the fictitious 
current q(Jp — Jn/b) is equal to the actual current times a linear func- 
tion of X. This function is always positive and varies from a minimum of 
1/6 to a maximum of 1. 
Combining (6.6) with (6.1) and integrating gives the equation 

that we had before. Now, however, E is not a minimum at the same point 
where s is a maximum. As before, when recombination is negligible 
throughout all of the intrinsic region, A determines the voltage; and, 
when recombination is important over part of the region, A determines 
both the voltage and the length of the intrinsic region Xp — Xa > 2L ^ 

i/gg- 

Combining (6.7) with (6.2) gives 

which is similar to the pre\nous (3.6) except that / is multiplied by the 
factor 1 + ^x/L, which varies from 1 + 1/6 to 1 + 6. The same argu- 
ments used in Section V apply here and show that the second term in 
brackets (the diffusion term) can be neglected except near the junctions. 
In other words, although / is always part drift and part diffusion, 
/(I + ^x/L) is approximately pure drift except at the junctions. 

Eliminating s between (6.9) and (6.8) and neglecting the diffusion 
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term in (6.9) gives the cubic equation 

^' _^('l + ^') = ('„-^Y-A (6.10) 



Ei^ aE\ LJ \2Li 

for the field distribution. 

In germanium, where b = 2.1, /3 = 0.406, Xp = 1.35L and x„ = 
-0.65Z/. The coeflicient of I/aEi therefore varies from 1.47 to 3.10, or 
by a factor of a little more than 2. This will introduce some asymmetry 
into the E versus x curve in the low field region where the fictitious car- 
rier flow Jp + Jn/b is by diffusion. It is evident that, as the voltage in- 
creases, the field versus x curve becomes increasingly symmetrical about 
the a; = point; so the effect of having b 9^ 1 is simply to shift the 
field distribution along the x axis. 

Recombination 

The arguments of section 4 again apply. Where recombination is im- 
portant, n — pis small compared to n -b p, so g — r = g{l — s^). The 
diffusion term dominates in the fictitious particle flow Jp + J,i/b; that 
is, E^IE\ is small compared to s, so (6.1) becomes 



-- = -2?iiZ) — 
ax 



J,+'-^= -2niD^ 



The continuity equations give 



So again we have 



dx' 2W 



(6.11) 



The solution joins the no recombmation solution where s = A — 
{x/2Lif. Therefore A is again related to So , the maximum s, by A = 
so(l - soVS) and the s versus x curve is given by (4.8) and is symmetrical 
about the point where s is a maximum. When the recombination solu- 
tion joins onto no-recombination solutions, there will be a different 
no-recombination solution on each side of the recombination region. 
The junctions will be at the points .t^ and .t„ on the respective no-recom- 
bination solutions. The length of the intrinsic region will not be Xp — 
x„ = 2L since the x = points are different on the two no-recombination 
solutions and are separated by a region of maximum recombination. 
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To find E when s is known we express the current I = q(Jp — J„) 
in tei-ms of s and E. Since n ~ p is small compared to n -\- p, we set 
n = p = sTii in. (2.2) and obtain 

Thus the current contains both a drift and a diffusion term. This is to be 
expected for unequal mobilities. When holes and electrons have the 
same concentration gradient, the electrons, which have the higher dif- 
fusion constant, diffuse faster than the holes; hence the diffusion gives a 
net current. It is seen that in the recombination region the total carrier 
concentration has a symmetrical distribution about the point where it 
is a maximum but the field remains unsymmetrical. 

Junction Solution 

When (Eo/EiY is large compared to I/aEi the junction solution is 
independent of 6; so the solution obtained in Section V is vahd. In all 
cases the junction solution can be found using the method of Section V. 
The effect of b will be small over most of the range where the junction 
solution holds because the concentration of one type of carrier will be 
neghgible. To be exact, / in (5.8) should be multiplied by the factor 
(1 -[- ^Xq/L), which can be taken to be (1 + h)/2h at the NI junction 
and (1 -f- &)/2 at the IP junction. Instead of equation (5.7) we have 

as can be seen by differentiating (6.10) with Ei = 2Z/,- = o- =1. 

VII. EFFECT OF FIXED CHARGE 

This section will deal briefly with the case where there is some fixed 
charge but where the carrier charge cannot be neglected. For no recom- 
bination, the field distribution is given by a first order differential equa- 
tion. Solutions in closed form are obtained for the case of pure drift flow. 
For recombination and charge neutrality the solution in Section IV is 
valid provided the fixed charge is small compared to 7i,- . We have seen 
that at large fields the E versus x curve becomes linear, correspond- 
ing to a fixed charge density of Ni where Ni = -v^2n,£/i,- . Thus 
the fixed charge may have a dominant effect on the space charge while 
ha\dng a negligible effect on the solution in the range where recombi- 
nation is important. 
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Let the density of fixed charge he N = Nd - N^ = excess density 
of donors over acceptors. A'' may be either positive or negative. In what 
follows we shall assume that A^ is positive. So the structure is NvP 
where v means weakly doped n-type. Equations (2.2) for the hole and 
electron currents remain unchanged. Poisson's equation become.s 

^= a{p-n-\-N) (7.1) 

dx 

We shall deal with the general case of arbitrary mobilities. As in Section 
VI it is convenient to deal with a fictitious current g{Jp — JJh) and a 
fictitious particle flow Jp + J^/h. The extra term in (7.1) drops out by 
differentiation when (7.1) is substituted into the equation for Jj, —J„/h 
so (6.2) remains unchanged. However, instead of (6.1) we have 

So the fictitious particle flow is no longer the gradient of a potential 
involving only E and s. 

No Recomhination 

As in Section VI the continuity equations can be inmiediately inte- 
grated to give (6.6) and (6.7). Again / is given by (6.9) where the dif- 
fusion term on the right can be neglected except at the junctions; so 
again we have asE = 1(1 + px/L). Substituting this into (7.2) and com- 
bining (7.2) and (6.6) gives a first order differential equation for E versus 
X. It is convenient to again use dimensionless quantities with Ei , 2Li 
and tjEi as the units of field, length and current respectively. Then the 
differential e(|uation becomes 



dx 
where 



r^;^ ^Ul + ^\j^ 2{x + aE) (7.3) 



N 

and as before Ni = V2ni£/Li , which is around 4 X 10'« in germanium 
at room temperature. The solution of (7.3) contains one arbitrary con- 
stant (which corresponds to A in the A'' = case). The lower hmit on 
the constant is determined by the necessity of joining onto a recombina- 
tion solution where s approached unity. The positions of x„ and .r^ of 
the Nv and vP junctions respectively are given by (6.5). 
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In the region of low fields where E is comparable to or less than /, 
(7.3) would have to bo solved graphically or on a machine. At higher 
fields the equation is easily integrated as discussed below. 

Case of Pure Drift 

When the flow is entirely by drift, E^ ^ I and (7.3) becomes 

dE X , ._ ,, 

^ = 2 + " (^■*) 

which is made integrable by the substitution E = yx. A family of solu- 
tions for positive E throughout the v region is 

{E - a.xy'iE + a^x)"' = ^o"'"""' (7.5) 

where 2ai = \/4 + a^ + a and 2a2 = -^74 -\- a' — a and Eo is the 
value of £■ at X- = 0. For an intrinsic region A'^ = a = and (7.5) reduces 
to E^ — Eq + x^, which is the same as (3.9) for negative A. Fig. 9 shows 
several curves for various values of -£"0 . These remain above, and at 
large distances approach, the asymptotic solutions E = UiX on the 
right of the origin and E = —azx on the left. These curvesdiffer from the 
corresponding curves for an intrinsic region in that the straight line 
asymptotes now have slopes of ai and —02 instead of ±1. Toward the P 
side the slope is greater because the positive change qN of the excess do- 
nors is added to the charge of holes. Toward the A'^ side of the v re- 
gion the slope is reduced because N compensates to some extent for the 
electron charge. As a increases and the v region becomes more n type, 
the solution approaches that for a simple NP jmictiou, where E = ax 
on the A'' side. 
Another set of solutions of (7.4) are given by 

(aix ~ ET'iatX + E)"' = ai'W'x,^ (7.G) 

Several of these are shown in Fig. 9. They remain below the linear 
asymptotes and go through zero field at a: = ±Xe . Actually these will 
join onto solutions of the more general equation (7.3) when E becomes 
small and the diffusion term becomes important. 

Recombinalion. When the fixed charge density is small compared to the 
intrinsic hole and electron density the treatment of recombination in 
Section IV remains valid. The recombination solution joins onto a solu- 
tion of (7.3) at small fields. When A'^ is comparable to », the recombina- 
tion solution is difficult even with the assumption of charge neutrality. 
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Fig. 9 — Field Distribution in the Range of Pure Drift for a fixed charge 
N = Ni , or a = 1. 
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APPENDIX A 



Prim's Zero-Current Approximation 

Prim's analysis is based on the assumption that the hole and electron 
currents are negligibly small differences between their drift and diffusion 
terms. Setting /p = /„ = then gives n and p as functions of the po- 
tential, which is found by substituting « and p into Poisson's equation 
and solving subject to the boundary conditions at the junctions. These 
conditions involve the applied bias and the majority carrier densities in 
the extrinsic regions. Since the current is assumed to vanish, the phe- 
nomena of carrier generation and recombination do not enter the 
problem and the results are independent of carrier mobility. The results 
will be exact when there is no applied voltage; the potential drop across 
the unit is then the built-in potential. In this appendix we use an internal 
consistency check to see for what values of applied bias the analysis 
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breaks down. First we find where the carrier concentration is in error by 

finding the bias at which the minimum drift current as calculated from 
qn(n + p)E becomes equal to the total current, as found from the excess 
of generation over recombination in the intrinsic region. We then go on 
to find where the error in carrier concentration gives a sufficient error in 
space charge to affect the calculation of electric field. As we shall see, the 
zero-current approximation gives too low a carrier concentration in the 
interior of the intrinsic region. This will lead to serious errors in the field 
distribution only if the space charge of the carriers is important. When 
the bias is sufficiently high or the intrinsic region sufficiently narrow 
that the intrinsic region is swept so clean that the carrier space charge is, 
in fact, negligible, it will not matter that the calculated carrier density 
is too low, even by orders of magnitude. In such cases, the electric field is 
constant throughout most of the intrinsic region. 

In the following we shall, for simplicity, take & = 1 and assume that 
the extrinsic regions are equally doped so that the problem is symmetri- 
cal. 

Carner Density 

We now find where, on the zero current assumption, the drift current 
becomes equal to the total current. This involves knowing only the 
carrier concentrations and the field -E, in the center of the intrinsic 
region, where the drift current gfi(n + p)Ei is a minimum. By symmetry 
11 and p are equal here and ?i = p = u; exp (~qVj2kT) where Va is 
the applied bias. The minimum field Ei is given by the total voltage 
drop V and the field penetration parameter ij, which is the ratio of the 
minimum field to the average field. Thus j? = 2LEi/V where 2L is the 
width of the intrinsic regions. The difference between V and Va is the 
built-in voltage (2kT/q)/(7i(N/?ii) where A'' is the m^ority carrier 
concentration in the extrinsic regions. We now have for the drift current 
in the center of the intrinsic region 

,.(„ + p)i?=5,i3g)jexp(-|^) (Al) 

We next find the total current from the excess of generation over re- 
combination in the intrinsic region. From the zero current assumption, 

np = n,-^ exp (~qVa/kT) is constant throughout the intrinsic region. 
Hence g — ris constant. So the current / = q(g — r)2L = qL(ni^ — np)/ 
Tti i is 



^-'-^[^ --(-©] ^^^ 



1282 THE BELL SYSTEM TECHNICAL JOUENAL, NOVEMBER 1956 

Equatmg this to the drift current (Al) in the center of the intrinsic 
region gives 



LV g7 , fqVa 



= ^.-^.csch(^l (A3) 



The error in carrier concentration is less for narrower intrinsic regions 
and lower biases. Thus (A3) gives a curve of L versus F„ such that the 
zero current solution gives a good approximation to carrier concentration 
for points in the VaL plane lying well below the curve. As expected, for 
zero bias, the soUition is good for any value of L. However, for a bias of 
several kT/g, the solution for carrier concentration breaks down unless 
Lisa, very small fraction of a diffusion length. 

Carrier Space Charge. 

In Prim's analysis the carrier space charge is so low throughout most 
of the intrinsic region that the field remains approximately constant 
and equal to Ei . However there must be enough carriers present that 
the drift currents of holes and electrons can remove the carriers as fast 
as they arc generated. In this section we ask where the space charge of 
the necessary carriers becomes large enough that its effect on the field 
can no longer be neglected. Let AE be the change in field due to the 
space charge in the intrinsic region (not counting the high field regions 
near the junctions). Unless AE is small compared to /?,- the neglect of 
carrier space charge will not be justified. We shall find the ratio of AE 
to Ei . 

If the field is to be approximately constant, then the hole and electron 
concentrations can easily be found from the hole and electrons currents. 
We shall deal with applied biases of at least a few kT/g, for which 
recombination is negligible and the total current / = qg2L = qniL/r. 
Since g — r = g is constant, the hole and electron currents are linear in 
X and, for constant field, are proportional to the hole and electron con- 
centrations respectively. Thus the net space charge of the moving 
carriers q(p — n) is proportional to x and varies from zero in the center 
of the intrinsic region to qy near the IP junction, where n is small 
compared to -p and the current flows by hole drift, so / = ?mP^.- • Thus 
the maximum charge is 1/y.Ei and the total positive charge of the car- 
riers on the P side of the center is IL/2tiEi . This gives a drop in field 

„ _ alL _ arii kT /^ 
~ 2^ ~ TWiLi" 
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Dividing by Ei = r]'V/2L gives 

^ = Jl_ (I^\ (A4) 

Setting A^" equal to some fraction, say 20 per cent of -E. gives a family 
of curves for V versus L, with ?/ as a parameter. Prim lias plotted such 
curves in Fig. 11 of his paper. His curves will be good approximations 
when V for a given L and j/ lies above the V given by (A4). 

Prim's results are expressed in terms of the parameters U = 
qV/2kT and L = 2L/£t where £e is the Debye length in the extrinsic 
material. £.e is given by the same formula as £ except that N replaces rii . 
Substituting these into (A4) and setting AE = Ei/5 gives 

L = 3.57— St/ (A5) 

n,£ 

Prim's U versus L curves will be accurate up to the point where they 
intersect the corresponding curves from (A5). For germanium a reason- 
able value of NLi/ni£ is about 10^. This says that Prim's curves go bad 
at about L = 10\ which would be about 2.1 X 10~^ cm in germanium 
at 300°C. 

Branching of the V versus L Curves 

An effect which does not emerge from the zero-current analysis is 
that V may have several values for the same L and ?]. In other words 
the V versus L curve for given tj will have more than one branch. Specifi- 
cally, there will be a single V versus L curve up to a certain L at which 
the curve splits into three branches that diverge as L increases. This 
may be seen as follows: Consider an intrinsic region that is long compared 
to the diffusion length. Suppose a bias is applied that is low enough not 
to appreciably affect the space charge and potential drop at the junc- 
tions. A current will flow and a proportional, ohmic voltage drop will be 
developed across the intrinsic region. If the intrinsic region is long 
enough, this ohmic voltage may become large compared to the built-in 
voltage before the voltage drop at the junctions has changed appre- 
ciably. In this range the field penetration parameter will be rising from 
zero to about unity as V increases from the built-in voltage and ap- 
pi'oaches the ohmic voltage. As the voltage continues to increase, the 
space charge begins to penetrate the intrinsic region and a majority of 
the voltage drop comes in the space charge regions. Let L be the ef- 
fective length of current generation. When L is larger than a diffusion 
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length but small compared to the length of the intrinsic region, then the 
voltage drop at the ends of the intrinsic region ^\■ill be proportional to 
L^ while the current, and consequently the minimum field, will be propor- 
tional to L. Thus 77 will be proportional to 1/L and will decrease as V 
increases and the region becomes more swept. Finally the two space 
charge regions meet; then tj rises again with V and approaches unity. 
Hence, for a given tj and length of intrinsic region, there will be three 
different values of V. For lower L the dip in the tj versus V curve will be 
less, and there will be only one V for some values of ??. Since ti starts 
from zero at the built-in voltage and approaches unity for infinite volt- 
age, there must be either one or three values of V for every i}. Thus 
when the V versus L curve (or in Prim's notation the XJ versus L curve) 
branches, it branches at once into three curves. Prim's plot gives the 
upper branch in cases where all three are present. 



